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Daniele of the NASA Lewis Research Center Dynamics and Controls Branch. 
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1.0 SUMMARY 


The objectives of this program were to adapt the General Electric 
Pynamlc Blade Row Compression Component Stability Model for use with NASA 
developed linearized stability analysis and frequency response analysis 
subroutines and to demonstrate that these analysis techniques, which are 
used primarily in control analyses, give results which are equivalent to 
the time-dependent solution of the nonlinear equations describing the flow 
in a compression component. 

Verification of the applicability of these techniques was carried out 
by comparing results with previously obtained NASA J8S-13 engine data and 
with previously conducted time-dependent studies using the above mentioned 
General Electric mode] . The equations of this DK>del were linearized about 
an equilibrium operating point resulcing in a matrix. A, known as the 
Jacobian matrix. Based on the first method of Liapunov, the stability of 
the system was determined by either examining the signs of the real parts 
of the roots of the A matrix characteristic equation using the alternate 
formulation Routh-Hurwltz stability criterion or by actually solving for 
the eigenvalues and examining the signs of the real parts of the roots. 

The frequency response aitalyses were conducted using the A matrix and an 
input matrix B in conjunction with the NASA supplied subroutines which 
determined the linearized system transfer function using Laplace transforms. 

The J85-13 clean- Inlet-flow surge line was determined for each of two 
engines: one at 80, 87, 9A, and 100 percent corrected speeds and the other 

at 87, 94, and 100 percent corrected speeds. The resulting 87 x 87 A 
matrices were tested for stability and the predicted surge line accurately 
matched the experimental data and the time-dependent model results except 
for the first engine at 94 percent speed. The cause of this anomaly is not 
understood. The stability characteristics of the other engine were accu- 
rately predicted. 

The linearized frequency response analyses were conducted on the first 
engine with inlet-total-pressure oscillations and with exlt-flow-ftmction 
oscillations over the 0.1 to 200 Hz frequency range for 25 frequencies. 

These results were comoared with time-dependent-solution results at fre- 
quencies of 30, 60, lOu, 150, and 200 Uz. It was found that essentially 
identical results were achieved when the time-dependent model total-pressure 
oscillations were equal to ±0.5 percent of the mean inlet total pressure or 
when the exit* flow-function oscillations were eoual to ±1.0 percent of the 
mean exit flow function. 

Stability predictions were made for nne-per-rev inlet total-pressure 
distortions at 80, 87, and 100 percent corrected speeds. A two-sector 
parallel compressor model was used and resulted in a 187 x 187 A matrix. 

At each speed, surge point predictions were made for total-pressure distor- 
tions produced by a 42 percent blockage screen with angular extents of 180®, 
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90*. 60*. 30* and by a 49 percent blockage screen with angular extents 

of 190*, 90*, onH 60*. The 180* angular extent stablllix analysis results 
were cohered to existing tine-dependent nodel solutions aad this conpari- 
son showed that e<|uivalent results were obtained. The stability analyses, 
in general, produced results which were in k- eplng with other parallel 
conpressor results. That is, the loss in surge pressure ratio Is over- 
predicted especially at low speeds (80 and 87 percent) where the speed 
lines have lower slopes. In an effort to inprove the distorted-flow sta- 
bility predictions, inter-stage axlal-ga:» crossflow fron the high total- 
pressure sector to the low total-pressure sector was sisulated at one axial 
location in the conpressor. The amount of crossflow was detewined by 
obtaining the best natch of the predicted distortion attenuation charac- 
teristics with the neasured distortion attenuation characteristics for a 
180* total-pressure distortion at corrected speeds of 80- and 100-percent. 
Crossflows for other extent discortlcHis were scaled fron these results 
based upon area ratios. It was found chat crossflows at rotor 1 exit of 
1.2 percent and 1.0 per«.r*tt of sector inlet physical flow at 80 and 100 
percent corrected speeds, respectively, gave the best natch to the neasured 
attenuation characteristics. Even though this nisnatching of stages niould 
significantly nove the operating point if the Irv total-pressure sector in 
nany cases, the nean conpressor behavior at surge (loss in surge pressure 
ratio) was not isproved and in sone cases was worsoied due to unrealistic 
sector speed line shapes. Use of a stage prediction nethod in such cases 
is necessary since the distorted sector nay never reac_ the clean-inlet- 
fiow surge line due co the stage nisnatching. The linearized stability 
and the frequency response analyses enployed in this study gave results 
which were equivalent to tine-dependent solutions of the nonlinear set of 
equations. Such analyses offer significant econonic advantage when linearized 
analyses are applicable. 
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2.0 IHTBODUCnOH 


One of the most significant paraaeters of concern In the design and 
operation of gas turbine engines is surge nargln. Aircraft engine designers 
are continually striving to provide adequate surge aaigln to insure success- 
ful cteady- and translqit-CTglne operation. Experience has revealed that 
surge uargln can be degraded by a variety of factors such as deterioration, 
control tolerances. Reynolds maber, stator vlsriggiiig, accelerations, 
inlet distortion, etc. Of these factors, one of the largest and nost 
variable cwnsaaers of surge aargin Is inlet distortion. 

Since it is not feasible to ezperiaentally detemine the response of a 
coapresslon coaponent or an engine to all possible coablnatloos of distortion, 
the need has evolved idierd>y the stability (surge) of a coapresslon coaponent 
can be analytically predicted. There have been a nnaber of recent studies 
using tlae-dep«idait aethods. References 1 and 2 report the results obtained 
using stage by stage aodels and Reference 3 reports the results obtained 
using a blade-row by blade-row aodel. As the aodels have be c o aa sore 
coaplex. it has been necessary to turn froa analog or hybrid coaputatlon to 
strictly digital coaputatlon. Even for very efficient solutions of the 
tiae-dependent equations such as that used in the Reference 3 study, coapu- 
tation costs represent a significant proportion of the effort. An approach 
(References 4 and 5) has be«i developed at NASA Lmiis Research Center %diich 
potentially offers significant cost reductions. 

This approach Involves linearising the nonlinear coupled systea of 
equations which describe the conpresslon STStes about an equilibrlua operat- 
ing point. This operation resulta in a natrlx. the characteristic equaticm 
of which can be examined for the sign behavior of the real part of its 
roots. Instability will occur lAen any root has a positive real part. 

Further, it has been shown that these instabilities occur only at the 
experimentally determined surge line. Based on the first method of Liapunov, 
the n' nllnear system of equations will be stable if the linear system of 
equations is stable. As explained in Reference 4. the stage by stage 
analytical representation of a compressor as given In Reference 1 was used 
and the resulting system of equations were then linearized. The stability 
of the resulting linear system which repr^ented a J8S-13 compressor was 
determined. The results showed good agreemoit between the analytically and 
experimentally determined surge lines. 

The objectives of the program v&re to adapt the NASA developed technique 
for predicting surge in a compression component for use with the exisiting 
General Electric Dynamic Digital Blade Row CcMapressicn Component Stability 
Model (Reference 3) and to perform an evaluation of the capabilities of 
this stability analysis technique by comparing these results to time- 
dependent solution results and test data reported upon in Reference 3. 

These objectives were to be met by determining the clean-inlet-flow surge 
points for a J85-13 turbojet engine known as the "Moss" engine at the 80. 

87. 94, and 100 percent corrected speed lines. Then the response of the 


3 



J8S-13 coopressioa conponent to Inlet-total-pressure and exit flow-function 
oscillations over the 0.1 to 200 Hz frequency range was to be deterwined. 
The aaplitude and phase angle results were to be coapared with the results 
of the Reference 3 tiae-dependent program at five frequencies in the 20 to 
200 Hz frequency range. The aethod was then to be extended to predicting 
the point of instability with distorted inflow conditions by using a systen 
of equations describing a two sector parallel conpressor. Stability pre- 
dictions for a ooderate level of distortion with angular extents of 30, 60, 
90, and 180 degrees and for a higher level of distortion with angular 
extents of 60, 90, and 180 degrees were to be conducted at 80, 87 and 100 
percent corrected speeds. Each of these cases was to be run with and 
without circirafvxential crossflow. 

In this report, the aethods used in the co>irse of the study are 
reviewed and the results obtained in the course of the stability and fre- 
quency response predictions are discussed. 
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3.0 ANALYTICAL METHODS AND J85-13 ENGINE COMPRESSOR MM)EL 


The analytical locmilation of the Dynaalc Digital Blade Row Coopression 
Component Stability Model and the methods employed in this study to deter- 
mine the stability of the system of equations are reviewed. Further the 
J85- 1 j engine which is being Pimulated, the source of the blade row charac- 
teristics, and the manner in which the compression system is broken into 
volumes are discussed. 


3.1 DESCRIPTION OF BLADE ROM MODEL 


This review of the Dynamic Digital Blade Row Compression Component 
Stability model is taken from References 3 and 6. The complete set of non- 
linear partial differential equations which describe the transfer and storage 
of mass, momentum, and energy within a fluid are called the equations of 
change. These equations have been Integrated once over an arbitrary vnlme 
of the flow system to .'btain the macroscopic balances for quasi ooe-diaoi- 
sional flow without heat transfer and are reproduced below in the form in 
which they are useu in the dynamic compression compooent model. 


3p 


k 


3t 





— (V - W 

^**i i+1) 


*o r Vfi 

Sc L 8o 8c 


i+1 '**i+l 




- Vl> ^ ^B 


\ l^'i^i - «i+l "i+1 + "pl 
k 


(1) 

( 2 ) 

(3) 


The subscripted variables on the right-hand side of the equations r<-.fer to 
quantities at the inlet (i) and exit (i+1) of the control volume. Variables 
on the left-hand side refer to volume averaged quantities, i.e., in gener- 
alized fora 


- ^ ^ dV 

^ / dV 


(A) 


The energy equation (Equation 3) was derived by combining the equation of 
change for energy and one of the thermodynamic T ds relationships. 
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This set of equations (other than being applicable to quasi one-dimen- 
sional flows without heat transfer and to a finite » but small volume) prop- 
erly and exactly describes the state of a fluid in motion. In order •'o solve 
Equations 1 through 3, it is necessary to supply the caloric and thermal 
equations of state and expressions for Fg (Blade Force), (Mean Pressure), 
and Sp (Entropy Production) . 

In the following paragraphs, the blade force, mean pressure, and entropy 
production terms are discussed. 

The blade force Fg of Equation 2 represents the blade force acting upon 
the fluid. The blade force can be determined through reference to the fol- 
lowing sketch: 



F 


B 


F tan S 
T 


DZ 


where 


(5) 


F 


T 



^i+l^i+l^^l-H 
r, + r 


lli!i"u 

i+l 


0 


( 6 ) 


and is derived from the Euler Turbine Equation. The direction of the lift 
vector is assumed to be 
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1^ 

2 


(8i + 


^i+1^ 


(7) 
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The drag force (Fq) Is obtained from the following equation 


F 


D 





( 8 ) 


which is based upon an analogy with the drag coefficient for duct flows. It 
is noted that the drag force acts in the direction of the inlet flow 
The term Fqz in Equation 5 is then obtained from the relation 

[cos(S. - f„)] 

F = F = (9) 

DZ D cos Boo 

where the numerator transforms the force (Fd) in the flow direction to a 
force along Che mean chord line perpendicular to the lift vector and the 
denominator transforms Chat resultcnt force to an axial force. 

It should be noted that in steady flow a raomentua balance, in goteral, 
will not give the same total-pressure rise per stage as does an energy bal- 
ance . The latter is assumed to be correct. The reason for this difference 
is that the proper direction of the blade lift vector is not exactly the 
arithmetic average of the flow angles (Equaticxi 7). Comparison of the 
steady-state momentum and energy balance solutions perait:> the deteminatimi 
of a small "correction angle" which can then be added to Sg, to give the 
proper lift direction. Hence, in actuality the lift direction angle is 
written as 


6« - 1/2 (G^ + (10) 

Equation 7, then, is replaced by Equation 10. 

The mean pressure of Equation 2 represents the average of the inte- 
gral of the pressure distribution over the lateral surface area of the volume 
element as sketched below: 
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Although an analytical expression for the mean pressure acting on lossless, 
blade free volumes In steady flow can be derived. It leads to redundancy In 
the system of equations describing unsteady flow. Therefore, based upon 
steady-state momentum-balance analyses, an approximate linear expression for 
calculating blade-free volume mean pressure as a function of area convergence 
and inlet and exit pressures his been established for the J85-13 compressor 
model. Its form in terms of a scale factor (FV) for zero-swirl free volumes 
is given by the following equation: 


° ^i ^ (Pi+1> / (1 + FV) (11) 

Similarly, an additional correlation was established for the non-zero 
swirl free volume between the IGV and the first rotor as a function of ICV 
exit air angle. It should be noted that the correlations are given in Ref- 
erence 3 and are unique to this formulation and the J85-13 compressor. They 
are probably valid only for other compressors with similar Mach numbers and 
area changes. 

For blade-row volumes, iiivestigations have revealed that a good approxi- 
mation for the mean pressure -s two-thirds the higher of the inlet or exit 
static pressures plus one third the lower pressure. Deviations from this 
approximation are accounted for in the lift direction correction angle. 

The term Sp in Equation 3 represents the total rate of irreversible 
conversion of mechanical to internal engergy and, in the case of this model, 
represents the entropy production due to blade row losses. It can be ob- 
tained from the expression: 


s . iR ' Ideal 

^ PC /?, 


^i). 


( 12 ) 
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where the ideal relative total-pressure ratio which accounts for the change 
in pitch line rad 
written as 
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Mj is equal to the ratio of the blade- row- ex it pitch- line wheel speed to 
the inlet relative stagnation velocity of sound (2ifNr£^j/axj) . In the case 
of a stator, the ideal relative total-pressure ratio is equal to one. ihe 
actual relative total -pressure ratio requires knowledge of the relative 
total-pressure loss coefficient which is defined as 


Pt 

Ti+1 




- Pf. 


ideal 
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(14) 


Equation 14 can be rewritten in the form 


pfi ; ■ pii / 
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- Od 


1 

1 

( 



y/(y - 1) 


(15) 


Hence, Equations 13 and 15 when substituted into Equation 12 provide com- 
plete definition of the entropy production term. 

As might be expected, the input to the program, in addition to physical 
speed, inlet conditions and compressor geometry, requires the relative 
total-pressure loss coefficient (u') for each blade row, the deviation 
angle (6) for the rotors, the correction angle (E^.) for both the rotors and 
stators. The relative total-pressure loss coefficient and the deviation 
angle are derived from stage-stacking results and are based upon clean- 
inlet-flow test data obtained by throttling at constant speed. The correc- 
tion angle is obtained by comparing steady-flow force and energy balance 
solutions on a blade-row basis. These parameters can be represented as 
functions of incidence angle and are input to the program in this manner. 
These parameters, in conjunction with the velocity triangles and other 
ancillary relations, permit the determination of the thermodynamics of the 
fluid at each station. Blade-free volumes are treated as lossless volumes 
with no imposed blade force; hence, the Fg and Sp terms of Equations 2 and 
3 are identically zero. 

Time dependent solution of the system of equations (Equations I through 
3 and the relations for Fg, and Sp) that comprise the dynamic digital 
comp’*ession component model is effected through a Taylor series which 
establishes the values of the three independent volume-averaged variables 
at the next increment in time. In the case of this model and ^th references 
to the left hand side of Equations 1, 2, and 3, the variables p, W, and 
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ps are ones for which a solution is sought. Solution is now straight- 
forward and will be illustrated for one variable - the voluae-averaged 
density. Considering that this method is applicable to any volume, the 
subscript "k", indicating the k-th volume will be dropped. The Taylor 
series for volume-averaged density correct to second order can be written 
as: 


p(t+At) = p(t) + 


3p(t) 

3t 


At + 


3^p(t) 


at 


(At)^ 

2 


(16) 


where: 

p(t) is established by t< nitial conditions or from the previous time 
step. 


^ from Equation 1 and differentiating this 


equation with respect to time yields: 


u- , . , / aw. aw.^, \ 

a p(t) _ ij X i+i ) 

2 “ V \ at ~ at ' 

(7 C 


(17) 


Examination of Equation 17 reveals that the right-hand side is composed 
of derivatives of station values of flow with respect to time. Since Equa- 
tion 2 will supply only the derivatives of the voliome-averaged flow with 
respect tc time, use of an interpolation scheme for obtaining station 
values from volume-averaged values will permit Equation 17 to be solved for 
the second partial derivative of volume-averaged density with respect to 
time. Equations 1 and 17 then can be substituted into Equation 16 to 
obtain the estimate of the volume averaged density correct to second order 
at the next incranent in time. Equation 17 implies that first derivatives 
with respect to time of a large number of terms (e.g., Fg, P^, and Sp) will 
be required. Although these expansions are lengthy, they can be derived in 
a straightforward manner and will not be reprodu^'ed here. Similarly , this 
technique can be used for the remaining two variables (W and ps) and can be 
continued from one time step to the next for the desired number of time 
steps. 


3,2 DESCRIPTION OF STABILITY ANALYSIS TECHNIQUE 

The following paragraphs explain the methods which are employed in the 
present study to determine che stability of the compression system model at 
various steady-state operating points on a speed line. First, let us 
rewrite Equations 1 and 2 using a slightly different notation and Equation 
3 by expanding the left-hand side and substituting Equation 1 as follows: 


10 



( 18 ) 


k - i <"i - “m) 

k 


Wk 


r '^i+i^Zi+i 


!o r v!i 

\ L 8o 


So 


' ^i+1 ^i+1 


■ ■ ^1+1^ 


] 


(19) 


^ ■ ”i+l®i+l ■ ®k”i ■" ^k'^i+l ] 


( 20 ) 


If we choose Pk, W^, and sk> (volume averaged properties) to be state var- 
iables, then the nonlinear Equations 18 through 20 can be written in the 
general form 


f(X) 


(21) 


where X is an n vector composed of state variables and €(X) is a continuously 
differentiable n vector. In the case of one general volume, the k-th. Equa- 
tion 21 would take the form of Equations 22 throu^ 24 ^ere the right hand 
sides of Equations 18 through 20 have been shown functionally in terms of the 
applicable state variables. 

^k “ ^1 ^\-2’^k-2’®k-2’\-l’\-l’®k-l’\’\’\’^fc+l’\+l’®k+l^ 

K " ^2 ^V2’\-2’ V2’^k-l’”k-l»\-l’\’\’®k’\+l’\+l’®k+l^ 


\ “ ^3 ^\-2’\-2’V2’\-l’\-l»®k-l’^k’\’®k’^k+l’\+l’®k+l^ 


It is noted that the above representation is applicab)e to either bladed or 
free volumes and may require state variable information from two upstream 
volumes and one downstream volume. The state variables from the adjacent 
upstream volume and the adjacent downstream volume are required because of 
the interpolation technique which is employed. The second upstream volume 
is required since stator exit conditions are functions of the rotor Incidence 
angle. In general, not all state variables as indicated by Equations 22 
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through 24 will be required. For a compression system represented by k vol- 
umes we can write the system of equatlcns In the following general form 
noting that the time derivatives will not be functions of all the system 
state variables. 


^1 



“1 


^2^^1’^1’®1*^2’^2*®2’ 

h 


fj (Pj^,Wj^,s^,P 2,W2*®2’ ^k’\’®k^ 

h 


• 

«2 


• 

^2 

= 

• 

• 


• 

• 


• 

• 


• 

• 


• 

• 


• 

K. 


• 

\ 


^3k-l^^l’”l’®l’‘^2’“2’®2 \’”k’^k^ 

K 


^3k^‘’l’”l’®l’^2’”2’®2’ ^k’\’®k^ 


Keeping in mind the notation of Equation 25, we return to the more simple 
notation of Equation 21. Now let us expand the nonlinear system of func- 
tions implied by Equation 21 about an equilibrium (steady-state) operating 
point on a speed line in a multivariable Taylor's series. A new variable 
Y=X-Xg, is introduced which shifts the origin to the equilibrium point. At 
this point the higher order terms are neglected and Equation 21 becomes 


Y = A Y 


(26) 
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where A is a 3k x 3k matrix known as the Jacobian matrix. The A matrix can 
be written as 


^^1 

'^1 


^^1 

^^1 

^^1 


^^1 

3f 

9Pj^ 

3W^ 

?Si 

9P2 

3W2 




3s, 



'^2 

3f2 

^^2 

'^2 





3Wi 

3Sj 

aP2 

3W2 

3s2 





^^3k ^f3k S^3k 

3Pj^ 3Wj^ 3Sj^ 


^^3k ^f3k 

3 ^ 


The formulation of the A matrix elements is described in Appendix A. 

The stability of the system of equations can be determined according to 
the first method of Liapunov (Reference 7) who showed that if the eigenvalues 
of the A matrix have nonzero real parts, then the stability of the nonlinear 
system of equations at the equilibrium point 2^ is the same as that of the 
linear system of equations at Y 0. Therefore, if all the eigenvalues of 
the characteristic equation of the A matrix have negative real parts, then 
the equilibrium point ^ is asymptotically stable. 

Since only the sign of the real part of the roots (eigenvalue) is impor- 
tant to the stability of the system of equations, significant reduction in 
computation time can be gained by determining only the sign of the roots 
rather than determining the values of the roots themselves. The stability 
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and frequency response analysis techniques anployed during the clean-inlet- 
flow studies, where an 87 x b7 A matrix was being examined, were based on 
the work of Seidel (Reference 8). The cha-acteristic equation of the A 
matrix is obtained by the method of Danilevsky (Reference 9) which reduces 
the matrix to Frobenius form. Then the signs of the roots of the charac- 
teristic equation are found by the alternate Routh-Hurwitz formulation 
method (Reference 10). Although the Kouth, Hurwitz, and the alternate 
Routh-Hurwitz formulation methods for determining the signs of the roots 
are all equivalent, the details of the calculational procedures are differ- 
ent and offer varying degrees of computational simplicity. However, accord- 
ing to Routh, stabill*'y is insured if all the coefficients of the charac- 
teristic equation are present and positive (necessary condition) and that 
all the terms in the first column of the Routh array (Reference 7) are 
positive (sufficient condition). 


Much technical effort was expended during the early stages of this 
program in developing techniques to avoid u«. '-utatlonal underflow and 
overflow problems which were experienced whet, nalyzing the 87 ^ 87 Jacobian 
matrix for system stability. Two precautionary measures were undertaken to 
avoid these numerical problmns which inevitably occur for models with 
numerous volumes. A unit length scale factor was imposed to optimize the 
range of matrix elements by keeping the elements as small as possible. 

This scale factor scaled all units of length in the A matrix and resulted 
in a span of values from 10® through lO-* (versus 10~^ through 10^^ prior to 
scaling). An overall scale factor was then applied to the matrix, posi- 
tioning the range as equally as possible about unity, i.e., (10~3 - lO^). 
This further reduces the chances of computational problems by numerically 
balancing the matrix. 


Problems were encountered with the Reference 8 methods during the 
distorted-inlet-f low studies where a 153 x 153 A matrix was involved even 
though double precision arithmetic was employed. Considerable light was 
shed on the problems by Wilkinson (Reference 11). "In general, the reduc- 
tion to Frobeius form (or the computation of the characteristic equation) 
is much le5 satisfactory than the reduction to tri-diagonal form. For 
many quite Harmless looking distributions of eigenvalues the Frobenius form 
is extremely ill-conditioned, and it is common for double precision arith- 
metic to be quite inadequate." Upon investigation it was found that the 
Aircraft Engine Group of General Electric had purchased a subroutine pac- 
kage* that would circumvent these problems and would only require sinsle 
precision arithmetic since, in this case, the roots are found directly by 
reducing the A matrix to Hessenberg form (Reference 11) . This subroutine 
package was used throughout the distorted-inlet-flow studies. 


Interpretational , Mathematical and Statistical Library-IMSL 
Interpretational , Mathematical and Statistical Libraries, In<_. 
Bellalre 

Houston, Texas 77036 
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3.3 DESCRIPTION OF FREQUENCY RESPCT^SE ANALYSIS TECHNIQUE 


The frequency response analysis technique of Reference 8 was employed In 
this study and is essentially a Laplace transfer function analysis of a set 
of linearized equations. 

Y - A Y + B U (28) 


R = C Y (29) 

These equations describe 1) the response of the state variables at the equi- 
librium operating point to input forcing functions and 2) tue output response 
in terms of a linear combination of the o.her variables. 

The homogenecus linearized system of equations describing the compressor 
determines the 3k 3k matrix and is described in Section 3.2 (Equacion 26). 
The state variable fector, Y is a ?h x i vector. U is a n >« 1 vector input 
forcing function and ^ is a 3k x n matrix of the form: 


^'l 

^^1 


9u, 


^u 

1 

2 

n 

“^2 


sf-, 




1 


n 


(30) 


3f 


3k 


9u, 


3f 

Tu 


3k 


The partial derivatives used to form the ^ matrix are given in Appendix B. 

R is s m 1 output vector and C is a m x 3k matrix expressing the linear 
relationship between the state variables and the output variables. 

The laplace Transform equivalents of the system equations represented by 
Equations 28 and 29 are 


s Y(s) = A Y(s) + B U(s) 

(31) 

R(s) = C Y(s) 

(32) 
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where s is the Laplace variable. Equations 31 and 32 can be written in 
transfer function form as 


R(s) . 

^ - C (s I - A)'^ B 


(33) 


where is the identity matrix. As it stands. Equation 33 requires a 
matrix inversion with a free parameter and is not compatible with digital 
computer calculations. The solution is readily ob«-ained, however, for any 
one or all elements of the matrix transfer function. 

For example the transfer function between the k-th element of R and 
the 1-th element of U is 


KR},. det (W) 

- k ^ - 

i(U}j “ det (s I - A) 


where W is an appropriate system Vronskian matrix. Details are given in 
Appendix C. 

Upon expansion, the numerator and denominator of the transfer function 
represented by Equation 32 are each polynomials in s. Using one of the 
subroutines in Reference 8, the transfer function frequency response is 
numerically evaluated by substituting s ju where i» is a given frequency 
and j is the iuaginary ntaaber /-I. 


3.4 LtSCRIPTION OF J85-13 COMPRESSOR MODEL 


The previously discussed stability analysis and frequency response 
techniques were used to study the characteristics of the J8S-13 turbojet 
engine. The compressor of a J85-13 engine has eight stages with a variable 
camber IGV and variable third, fourth, and fifth stage bleeds located in 
ea<.h stator channel at the casing wall. Ihe IGV trailing edge flaps and 
the bleeds are ganged together and are scheduled as a function of corrected 
speed biased by compressor-face total temperature. A cross-section view of 
the J85-1? engine in the compressor and combustor regions is shown in 
Figure 1. The nominal IGV and bleed schedules for the bill of material 
engine are given in Figure 2. 

It is appropriate to discuss the two engines modeled in this study. 

The pure total-pressure distort: jn patterns that were simulated during this 
study were obtained during testing of a J85-13 engine known as the "Moss" 
engine (Reference 12). This engine was run in support of the NASA casing 
treatment program. The clean inlet and distortion data utilized in this 
model were obtained from the untreated onfiguration with solid compressor- 
case Inserts. It should be noted that this engine tends to be representative 
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Figure 1. Section View of J85-13 Compreeeor 




of a deteriorated engine. For conparative purposes, clean inlet flo«f 
studies were also aade on a J8S-13 engine known as tl>e "Mehalic" engine 
(Reference 13). NASA had used this engine to obtain pure total-tenperature 
distortion pattern data and the coobined total-pressure and total-teaperature 
distortion pattern data. 

For the purposes of this study, the coapressor aodel includes volumes 
upstream of the IGV to the distortion measurement plane and downstream of 
the OGV continuing to a choke plane located at the turbine nozzle diaphragm 
(A 4 ). The purpose of including these eztra volwes is to insure that 
realistic boundary conditions can be imposed. 

The compressor model consists of twenty-nine volumes. There are 18 
bladed volumes (one blade row per volume) consisting of the IGV, rotors 
1-8, stators 1-8, and the OGV and 11 free volimes. These free volumes 
consist of two volumes between the instrumentation plane and the leading 
edge of the IGV, a volume between the trailing edge of the IGV and the 
leading edge of rotor 1 , and eight volunes between the trailing edge of the 
OGV and the turbine diaphragm. This configuration is shown schematically 
in Figure 3. It should be noted that the length of a rotor blade row 
extends from the trailing edge of the upstream stator to the leading edge 
of the downstream stator and includes the axial inter-blade row gape while 
the length of a stator blade row extends from the leadii^ edge of the 
stator to the trailing edge of the stator. For further geometrical details, 
the reader is referred to Reference 3. 

The blade work on the fluid is accomplished in a distributed, but un- 
specified, manner across a rotor volime length. All losses are assumed to 
take place in rotors, that la, no losses are accounted for in blade free 
volumes or stator blade voliaes. Bowever, this is not a restriction of the 
model. The rotor deviation angles vary as functions of the incidence 
angle. 

At this point, it is worth mentioning that because a blade row fonni- 
lation is being used, stage characteristics Information was used in the 
form of a relative total-pressure loss coefficient and a deviation angle 
rather than the more often used non-dimensional work and pressure coeffi- 
cients as a function of flow coefficient. This approach tends to decouple 
the inlet and exit stations of a blade row in a dynamic analysis since 
volume storage of mass, momentw, and energy are permitted within the work 
producing volume. Further, splitting a stage volume into two blade row 
volumes will potentially double the frequency capabilities of a dynamic 
compression component model. The relative total-pressure loss coefficients 
and the deviation angles that were used during this study are given in both 
graphical and polynomial form in lleference 3. 

The third, fourth, and fifth stage bleed flows in the BK>del are 
removed at the exit of the stator while holding exit air angle constant. 

The percentage of inlet flow that is removed from each stage is given in 
Table 1. 
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Table 1. Bleed Reaoval Schedule (Percent 
of Inlet Physical Flow). 


Corrected 


Engine 

Speed 

Stage 3 

Stage 4 

Stage 

Moss 

80 

3.95X 

4.90Z 

5.78X 

Moss 

87 

2.28 

2.83 

3.33 

Moss 

94 

0 

0 

0 

Moss 

100 

0 

0 

0 

Mehalic 

87 

1.20 

1.49 

1.76 

Mehalic 

94 

0 

0 

0 

Mehalic 

100 

0 

0 

0 


During NASA testing of the Hoss and Hehalic engines, the Tx 2 control 
bias as shovn in Figure 2 was subverted. The Moss engine was controlled to 
294.3* K (70* F) or higher tenperature curve and the Mehalic engine had a 
Bodlfied schedule. Further details are given in Reference 3. 

The boundary conditions laposed upon the aodel consist of specifying the 
total pressure and total teaperature entering volisne 1 and specifying the 
value of the exit flow ftmction (Vi^g/Pg) at the exit of the last voluae. 
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4.0 RESULTS 


The results of the clean-ialet-flow, frequency-response, and distorted 
inlet-flow analyses of the J85-13 compression system are presented in this 
section. The clean-inlet-flow stability analyses vere conducted using the 
Routh-Hurwitz stability criterion, that is, only the signs of the roots of 
the characteristic equation were examined. The frequency response of the 
compressor was determined by evaluating the Laplace transfer function of 
the linearized system of equations. The distorted-inlet-flow stability 
analyses were accomplished by first constructing a parallel compressor 
model and then conducting the stability analyses by actually determining 
the roots (eigenvalues) of the characteristic equation. The signs of the 
real parts of the roots were examined (all negative values for stability) 
to determine if the operating condition was stable. 


4.1 CLEAN-INLET-FLOW SURGE LINE ANALYSES 


The first step in verifying the state equations in their linearized 
form was accomplished by duplicating the J8S-13 clean- inJ st-f low compressor 
map. Figure 4 provides a detailed view of the stability analysis process 
for the 80Z corrected speed line. These linearized analyses (L-A) results 
are compared to the results of the time-dependent (T-D) dynamic analysis 
program (Reference 3) and to engine test data (Reference 13). The Routh- 
Hurwitz stability analysis technique was used to verify the stability for 
five surge-free conditions on the speed line and to predict the surge point 
within a weight flow uncertainty of ±0.05 kg/sec (±0.11 Ibm/sec). Figure 5 
shows the clean-inlet-flow compressor map for the J85-13 "Moss" engine. 

The 80, 87, 94, and 100 percent corrected speeds are shown extending over a 
range of pressure ratios from surge to the region of the nominal sea-level- 
static operating line. 

At 94 percent corrected speed, the linearized analysis failed to match 
the surge point. At this point it was appropriate to question whether this 
failure is a result of a limitation of the linearized blade row model or is 
the result of some minor error in the Jacobian matrix formulation which had 
not been discovered. In an effort to determine the reasons for this failure, 
two areas of concern were investigated: 1) the formulation of the A matrix 
elements and 2) the blade characteristics input (loss coefficients, deviation 
angles, and lift direction correction angles) for the Moss engine 94 percent 
corrected speed line. The A matrix elements were independently rederived 
and to our best knowledge are error free. Further, as explained in the 
next section (Section 4.2) it was possible to obtain an independent assess- 
ment of the accuracy of the formulation by comparing the frequency responses 
of the state variables p, w, and s at selected stations with the time- 
dependent solutions. This assessment showed that close agreement could be 
obtained between the linearized analysis and time-dependent-solution state 
variables. In the process of making minor changes or revisions to the A 
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Figure 3. "Mohs" Engine Stability Analyses Results. 



matrix, it was noticed that the location of the 94 percent corrected speed 
instability point would move significantly although the location of the 
stability points of the 80, 87, and 100 percent corrected speed lines moved 
almost imperceptibly. At this juncture it was questioned whether the blade 
row characteristics might play a part in this sensitivity. The pressure 
coefficients 'i' (in part from which the blade row characteristics are 
derived) for each stage of the Moss engine were examined as a function of 
flow coefficient. As can be seen from examination of Figure 6 (taken from 
Reference 3), Stages 2, 3, and 4 are operating on the positive slope por- 
tion of their characteristic at instability (lowest shown flow coefficient). 
Another J85-13 engine known as the Mehalic engine had been modeled at 94 
percent corrected speed during the Reference 3 study. Examination of the 
pressure coefficient - flow coefficient relationships for this compressor 
as depicted in Figure 7 reveals that only two stages are operating on the 
positive slope portion of their characteristics at instability. It should 
be noted that in general, the fewer the number of stages operating on 
positive slope characteristics, the more stable the compressor. The N^alic 
engine clean- inlet.-f low characteristics were analyzed using the same line- 
arized technique as applied to the Moss engine. The results are shown in 
Figure 8 and indie..'. ' good agreanent with test data and the time-dependent 
model. One can question as to whether the "Moss" engine 94 percent corrected 
speed stage characteristics are realistic representations of the engine. 

Since these characteristics were used in the time-dependent model to obtain 
an accurate prediction of the surge line, it is assumed that there is 
nothing inherently wrong with these characteristics. Thus, there appears 
to be no basic error in the formulation of the linearised model, but there 
does appear to be a sensitivity to the number of stages operating on the 
positive slope of their characteristics. 

These lines of investigation have provided Insight, but not the reason 
for the linearized analysis to Inaccurately predict the point of instability 
for the Moss engine at 94 percent corrected speed. It is suspected that 
the reason for this discrepancy involves two factors: 1) some of the 

assumptions anployed to permit linearizing the time-dependent model or some 
minor undiscovered error in an A matrix element (s) act to subtly decouple 
the stages and 2) the number of stages which operate on the positive slope 
portion of their characteristics. Thus, these two factors acting in con- 
cert tend to create a situation which is less stable than should be antici- 
pated. 

During the checkout phase when it was necessary to verify the accuracy 
of the derivation of the A matrix elements, it was found that the lack of 
output information which could be compared with an Independent solution 
seriously impeded the debugging process. Clean-lnlet-flow stability checks 
in essence provided only "go" or "no go" type of information and did not 
provide detailed dynamic responses of thermodynamic or flow parameters 
(state variables) which could be compared to the Reference 3 time-dependent 
model. However, the frequency response analyses discussed in Section 4.2 
provided state variable information at each station within the model. By 
selectively exanlnlng the state variable responses in blade free volumes. 
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h'iijure 8. "MehaJic" Engine Stability An.ilyses Results. 



rotor blade row volinnes, and stator blade row volumes, and comparing these 
responses with time-dependent solution results, it was possible to debug the 
A and E matrix element formulations in a more efficient and orderly fashion. 

For surge line determination alone, the linearized stability analysis 
method costs are approximately one-flf teenth of those incurred by the time- 
dependent method of Reference 3. 


4.2 FREQUENCY RESPONSE ANALYSIS 

The frequency response? of the "Moss" J85-13 compressor wete generated 
using the Jacobian (A) and input response (B) matrices and 'he Reference 8 
digital subroutines to determine the Laplace transfer function. The fre- 
quency response of compressor-exit state variables (p, W, s) to variations 
in compressor inlet total pressure and exit flow function over a range of 
0.1 to 200 Hz was determined. All investigations were carried out at the 
100% corrected speed operating point [PR=6.636, Wc=19.68 kg/sec (93.3 
Ibm/sec)] for 25 frequencies. The time-dependent program (Reference 3) was 
run concurr-'^ntly to provide comparative cases at frequencies of 30, 60, IOC, 

] 50, and 2GJ Hz. Combined plots of amplitude ratio (output to input) and 
'-base angle are given for the range of frequencies Investigated. Figures 9 
through 11 show the compressor exit state-variable frequency responses to 
inlet total-pressure oscillations. Good agreement In both amplitude and 
phase angle is noted between the linearized and dynamic program rest Its. 

The time-dependent program total-pressure oscillations were ±0.5 percent 
about the specified inlet pressure. It was found that time-dependent program 
oscillations of ±3.0 or ±2.0 percent of the Inlet total prest»_ were too 
large to obtain a match between the time-dependent program and the linearized 
analyses. Hence, a compression component Introduces significant nonllnearlties 
at low amplitudes. Application of the linearized technique for determining 
planar wave dynamic distortion transfer (Reference 6) must be done with 
great care since an engine In a highly manueverable aircraft may be subjected 
to inlet-produced unstec.ly flows where the ampli.tudes may be on the order of 
10 percent of the inlet mean total pressure. 

The amplitude and frequency responses to exit- flow- function oscillations 
are shown in Figures 12 through 14. The time-dependent program exit-flow- 
function oscillations were maintained at ±1.0 percent about the mean exlt- 
f low-function value. With this type of boundary condition, it was found 
that che dynamic program settled out slowly from its "start up" transient. 
Because it was economically impractical to let every dynamic case completely 
settle out, some of the minor differences between the linearized and dynamic 
predictions of the amplitude and frequency responses of the lower frequencies 
may be attributed to the fact that a stationary state had not been obtained. 

A detailed analysis of both the Inlet total-pressure and exlt-f.low- 
functlon oscillation cases shows that as 200 Hz is approached a systematic 
divergence develops between linearized and time-dependent program analyses 
amplitudes. At this time, no explanation has been found which satisfactorily 
explains this divergence. It may be that at frequencies equal to or greater 
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Figure 10. Compressor- Exit-Plow Frequency Response to Inlet 'Total-Pressure Oscillation 




Figure 12. Campreeeor«Exit-Oenelty Frequency Reeponee to CombiiMtor 
Exit Flow»Functlon Oscillation. 
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than 200 Hz for the J85-13 coopression systen, even though only very snail 
perturbations are being inposed upon the tiae-dependent nodel, the nonlinear 
terns in the atacrobalances (Equations 1 through 3) introduce a noticeable 
contribution to the state variables. 

It was concluded that the linearized frequency response and the tine- 
dependent BKKlel give essentially equivalent results (in terns of the state 
variables) in the frequency range of najor interest. It is now possible to 
conbine the state variables at the coaq>ressor exit to obtain the response of 
a nore easily recognized and neasurable quantity, the total pressure. This 
is accostplished by deriving the functional dependence of the total-pressure 
fluctuations to the fluctuations of the local state variables. Total deriva- 
tives are used to obtain this relationship. Static pressure can be expressed 
in terns of static density and entropy as: 


P 


/P_|^ ^(s-Sj.)(y-1)/R 


(35) 


where the subscript "r" denotes reference conditions. 

Sinilarly, Mach nunber can be expressed in terms of static pressure, 
static density, and weight flow as 


m2 


W2 

PA^ igo P 


(36) 


Total pressure can be related to static pressure and Mach nunber at a point 
by the isentropic relationship 


P I 1 + ^ m 2 j ^ ^ 


The total derivatives of Equations 35, 36, and 37 are respectively: 


dP = P 


(1 ds) 


(37) 


(38) 
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= 2WpPdW - w2 (pdP 4 Pdp) 
p2 a2 

dP? = (l + ^ m 2|^"^ dP + p(^j |i + m 2 ^ dM^ 


(39) 


(40) 



Substituting Equations 38 and 39 into Equation 40 yields the total- 
pressure perturbation response characteristics about a steady-state operating 
point in terns of local state-variable perturbations. In this case, the 
differentials dp, dW, and ds represent the frequency response (anplitude 
and phase) of these state variables in coaples fora. It is now possible to 
explicitly calculate the response of conpressor-exit total pressure and 
directly conpare the results to those obtained froa the tine-dependent 
program. Figures IS and 16 show this comparison at the lOOZ corrected 
speed operating point for inlet total-pressure and exit-f low-function 
oscillations. Shown in the figures are normalized pressure responses to 
input nomalized inlet total pressure and input normalized exit ilow func- 
tion, respectively. As expected, the compressor attenuates the inlet 
total-pressure fluctuations. In conclusion, it has been shown that the 
linearized J85— 13 frequency response characteristics match the results 
obtained froa the Dynamic Digital Blade Row Compression Component Stability 
Model when small amplitude boundary conditions are imposed. 

However, it is apparent that for even mr derately low amplitudes, 
the effects of the nonlinear terms as included in the time-dependent model, 
can cause significant deviation of the time-dependent results from the 
linearized analysis results. If one were interested in large amplitude 
planar-wave transfer such as encountered with inlet buss, then it would be 
necessary to obtain time-dependent solutions of the nonlinear equations 
rather than a linearized analysis. 

It is important to note that the linearized analysis frequency response 
technique costs are independent of frequency and in fact, the frequency 
responses at many frequencies can be obtained for essentially the same cost 
as one frequency. Hotter, the cost of obtaining frequency responses is 
directly proportional to the nueber of stations idiere the frequency response 
is desired. On the other hand, the dynamic program is frequency sensitive 
because of the time to settle out and the period of a given wave. For a 
given configuration, the costs are Ind^endent of the nisaber of frequency 
responses desired. 

It is clear that if the linearized frequency response of a compression 
component at selected stations is desired, than the linearized technique of 
this section offers considerable economic advantages. However, if the 
frequency response of each stage at a given frequency is required, then the 
dynamic BK>del offers economic advantages for frequencies greater than 
approximately 20 Hz. 


4.3 DISTORTED-INLET-FLOU ANALYSES 


The blade-row parallel compressor model of Reference 3 was put in 
linearized form and resulted in a 153 x 153 matrix for a two sector model. 
The model was modified to permit flow at a given stage to be removed from 
the high total-pressure sector at its total pressure and to be reintroduced 
to the low total-pressure sector at its total pressure assuming that the 
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loss In total pressure occurred as a result of the crossflow In the Inter- 
stage axial gaps. The results of the distorted flow analysis with and 
without crossflow are discussed in the following paragraphs. 


4.3.1 Parallel Coapressor without Crossflow 

During the NASA casing treatment program the J85-13 engine was tested 
with 1/rev circumferential total-pressure distortion generated by 42 percent 
(7-1/2 M) blockage screens with angular extents of 180**, 90°, 60°, and 30° 
and with 49 percent (9M) blockage screens with angular extents of 180°, 90°, 
and 60°. The results of this testing are reported in Reference 14. 

These data were simulated using a two-sector parallel compressor model. 
The stability of the equations was determined using the A matrix formulation 
and actually calculating the roots of the characteristic equation. The 
roots were then examined to determine if a region of instability had been 
encountered as characterized by the presence of a positive real part of a 
root. The stability test was carried out for two levels of distortion 
(Table 2) at 80, 87, and 100 percent corrected speeds for the "Moss" engine. 
The stability tests were not carried out at 94% corrected speed since unex- 
plainable results had been obtained at this speed line during clean-lnlet- 
flow stability predictions. 

The results of the nnderate level of distortion (42 p<?.rcent blockage) 
stability analyses are given in Figures 17 through 20 while the results for 
the high level (49 percent blockage) of distortion are given in Figures 21 
through 23. The 180° 1/rev circumferential total-pressure distortion 
results of Figures 17 and 21 can be compared directly with the time-dependent 
model analysis of Reference 3. These comparisons show that the two methods 
give equivalent results. Examination of Figures 17 through 23 shows that 
the stability analyses of the two-sector parallel compressor model gave 
results which accurately reproduced the experimental speed lines at 80 and 
87 percent corrected speeds, but at the higher level of distortion at 100 
percent corrected speed, the flow shirt to the low side was not being repro- 
duced. All the results of Figures 17 through 23 show the typical classical 
parallel compressor results as far as predicting loss in surge pressure 
ratio, that is, it is predicted most accurately at high speeds where the 
speed lines are steep and it is overpredicted at the lower speeds and for low 
extent distortions. 

The stability analysis method employ > x>. ;nxs parallel compressor 
study offers economic advantages over the time-dependent model since its 
operating costs are approximately one-fifth of those incurred by the operation 
of the time-dependent model. 


4.3.2 Parallel Compressor with Crossflow 

In an effort to overcome the deficiencies of the parallel-compressor 
approach, that is, the overprediction of the loss in surge pressure ratio, 
it was hypothesized that inter-blade-row axial-gap crossflows might be the 
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Table 2. Measured Total-Pressure Distortion Levels 


42% Blockage Screen 49% Blockage Screen 


% N//e 

Extent 

^FT/PT>„ax* 

APT/PT>Avg** 

APT/FT>Max 

^T/PT>A^g 

100 

180 

0.1118 

0.0559 

0.1367 

0.0683 


90 

0.0918 

0.0689 

0.1114 

0.0835 


60 

0.0870 

0.0725 

0.1063 

0.0885 


30 

0.1009 

0.0925 

— 

— 

87 . 

180 

0.0594 

0.0297 

0.0677 

0.0339 


90 

0.0476 

0.0357 

0.0614 

0.0461 


60 

0.0446 

0.0372 

0.0558 

0.0465 


30 

0.0511 

0.0468 

~ 

— 

80 

180 

0.0423 

0.0212 

0.0474 

0.0237 


90 

0.03S1 

0.0263 

0.0428 

0.0319 


60 

0.0331 

0.0276 

0.0388 

0.0323 


30 

0.0373 

0.0342 




^PT/PT>Max “ Max " ^T MinJ/^T Avg 
APx/P-j)^^g “ (Pt Avg - Pt Min^/^T Avg 


41 



L-A Stability Limit 



i2 


Figure 17. 180 , l/Rev Total-Pressure Distortion ( 42 % Blockage) Stability Analyses 




Figure iH. 90**, l/Rev To tal>Pressure Distortion (42% Blockage) Stability Aiialyse 






NASA Test Data 
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Corrected Plow, kg/sac 


NAbA Test l>uta 
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Figure 21. l80 , l/Hev Tola l-PreH«<uro HiMtortion (^*9% Dlocknge) Stability Analyses 



O NASA Test Data 
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Figux-e 22. 90 , l/Rev Total-Pr«*t<ur« Distortion (^% Blockage) stability Analyaea 
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■echanim responsible for the loss in surge pressure ratio being less than 
predicted. The anount of crossflow was selected based upon natching neasured 
intemal-stage-distortion ai^lif icacioos (Beference 15) at 80 and 100 
percent corrected speeds for the 180° 1/rev high total-pressure distortion 
as shown in Figure 24. It was found that a crossflow of 1.2 and 1.0 percent 
of the undistorted sector inlet physical flow between stator 1 and rotor 2 
at 80 and 100 percent corrected speeds, respectively, b^t watched the 
stage anplif ication data at the given operating point. The crossflow at 87 
percent corrected speed was assmed to fall between the 80 and 100 percent 
values uid was assiHKd to be 1.1 percent of sector inlet physical flow. 
Initially, the percentage crossflow for the other extents was established 
by nultiplying the percentage crossflow at 180° by the ratio of the angular 
extent being sinulated to 180° (ratio of the areas). However, this led to 
a speedline shape for the 30° extent distortion at ICO percent corrected 
speed which significantly departed fron the experinoitally detemined 
speedline shape. The percentage crossflow was lowered for the 30° extent 
distortion until the speedline regained a proper shape. The percent cross- 
flow at ocher speeds was lowered proportionately and the 180° and 30° 
points were connected by straight lines. This crossflow schedule as a 
percentage undistorted sector inlet physical flow is glvoi in Figure 25. 

The results of Che crossflow studies are shown in Figures 26 through 
29 for Che 42 percent blockage screois and in Figures 30 throngjh 32 for the 
49 percent blockage screens. It is Interesting to note that at 80 and 87 
percent corrected speeds with the 42 percent blockage screens, only an 
insignificant increase in Che predicted surge pressure ratio levels as 
coepared to the levels without crossflow were noted. At 100 percent cor- 
rected speed, the introduction of crossflow woreene die data natch and in 
fact causes an unrealistic hooking over of the wean speed line prior to 
indicating instability. This is due to the fact that stage one is forced 
to operate on the very low flow portiim of its characteristics. 

Bsseatially the sane c oion ts can b» jtated about die 49-percent- 
blockage-scrd data except that at the low corrected speeds, the pressure 
ratio at instability is slightly lower. 

In sunwary, it does not appear that the introduction of crossflow and 
its attendant nlswatching of stages resulted in a better prediction of the 
loss of surge pressure ratio due to total-pressure distortion. It is 
possible that this result is biased by assuning a constant percentage of 
crossflow indepidenc of operating line and/or by assuming that the shape 
of the attenuation characteristic (Figure 25) remains unchanged for extents 
of distortion less than 180 degrees. 


49 





NASA Test Data 
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Figure 2(>» l80 , 1/Rev Total-i^reHsure Distortion (U2% Blockage) With Interstage Crossflow 

Stability Analyses. 



o NASA Teat 
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Figure 27. 90°, l/Rev Total-Preeeure Distortion (43% Blockage) With Interstage Crossflow 

Stability Analyse'. 



ONASA Test Data 
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Figure 28, 6o°, l/Rev Total -Pressure Distortion ( 42 % Blockage) With Interstage Crossflow 

Stability Analyses. 
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Figure 29. 30°, l/Rev Total-Preasure Distortion (42?4 Blockage) With Interstage Crossflow 

Stability Analyses. 



O NASA Test Data 
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Fifjure 30. l8o , l/Kev Total-Pressurt ,-iStortion (49% Blockage) With Interstage Crossflow 
Stabilit> Analyses. 



L-A Stability Liait 
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Figure 31. 9o”, l/Rev Total-Pressure Distorti''* < ’ 'lockage* With Interstage Crossflow 

Stability Analyses. 



^ l-A Stability Limit 



Figure 6o , l/Rev Total-Preeeure DiMtortion Blockage) With Intoratage Crowaflow 

Stability Analyaea. 







5.0 OQMCmSIOIIS AMD BECOWgllDATIOMS 


This progrsB is an esaople of the cross-fertilisation of technology 
areas which saist take place if progress is to continue and if a greater 
understanding of physical processes is to be <d>tained. In this case* lin- 
arized stability and frequency response analysis nethods idiich are used 
prlaarily in control analyses are applied to the detemination of gas- tur- 
bine-engine conpressor-aerodynanlc stability and frequency-response charac- 
teristics. The following conclusions can be drawn fron the results of this 
study: 

o The linearized stability analysis based upon the first zetbod of 
Liapunov using the Jacc^ian A natriz and the alternate fonulatlon 
of the Routh-Hurwitz stability criterion accurately reproduces 
the experinentally deteinlned J85-13 clean-inlet-flow surge 
lines. These surge line results are also in close agrconent 
(±0.05 kg/sec) with the results Stained fron the tlM-dependent 
nodel of Reference 3. 

o At 94Z corrected speed for the Tfoss" eDglne* surge is predicted 
at a pressure ratio lower than detemined experiBentelly or 
predicted by the tisie-dependant K^el. Althon^ greet care was 
taken to insure accursi^* e SBell forsnlatlon error in tte A 
■atrix or a nlnor progrvBing error cannot be ruled out at this 
tine. Further, it is possible dmt nauspcloos introduced to 
pemit linearising the tine-dependent nodel act to subtly decouple 
the stages and this dacoupli^ is snpllfied by the dieracter of 
the stage characteristics. 

o The llnaarisad-stabillty-analysis nethod offers clear econonlc 
advantages over the thze-dependent nodel for clean-inlet-flow 
surge line deteminatlons since its associated costs are one- 
fifteenth of those associated with the tine-dependent nodel. 

o The linearized-frequaicy-analysis technique produced results 
which accurately matched the tine-d^^ndent nodel up to nearly 
200 Hz. 

o Since it was detemined that the linear analysis appeared to 

apply only for total-pressure oscillatious on the order of ±0.5 
percent of the mean inlet total pressure, great care mist be 
exercised if the nethod is used to detemlne planar-wave dynamic 
distortion transfer coefficients since inlets of highly maneuver- 
able aircraft nay produce amplitudes on the order of ±10.0 percent 
of the mean inlet total pressure. 

c Depending ipon the frequency of interest and the number of oLations 
where frequency rwponses are desired, there are economic tradeoffs 
which should be made before choosing either the linearized analysis 
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techolque or the tiae-dependent ■sthod. In general if one or two 
discrete frequencies above 20 Hz are of interest, the tlaa- 
dependent nethod is nore econoaical. However, if a coaplete 
frequency response is of interest, then the linearized frequency 
response analysis is nore econonical. 

o The linearised stability analysis of a two-sector blade row 
parallel cowpressor aodel gave results for 180* 1/rev inlet 
total-pressure distortions which reprodnced the results obtained 
froa the t la c -dependent aodel. 

o The linearized stability analysts of the two-sector blade row 

parallel coapresaor aodel for the 180*, 90*, 60*, and 30* dlstor- 
tioc extent:; give results which are typical of parallel coapres- 
sors, that is, the loss in surge pressure ratio is oven>redicted 
especially for the lower extents of distortion and for the lower 
spee^ where the speed lines are of lower slope. 

o Allowing for crossflow froa the high total-pressure sector to the 
low-total pressure sector Just downstreaa of stage 1 did not 
laprove the accuracy of the loss in surge pressure ratio. In 
soae cases the saall aaount of crossflow showed sufficient stage 
nisaatcb to force an unrealistic "hooking over" of the speed line 
prior to instability although stability predictions were aade in 
this region. 

o The costs of asking a linearized stability analysis predictions 
of a two-sector parallel coapressor are approxiaately one-fifth 
of those incurred \xf operating the tine-dependent aodel. However, 
as the niaber of sectors increases, the A aatrix used in waking 
the linearized stability analysis becjaas large very rapidly. 

Thus, when waking a multisector parallel cfMupressor analysis. It 
is necessary to evaluate idiether the tiae-d^>endent approach aay 
offer a nore econonical approach. 

The linearized analysis stability prediction and frequency response 
techniques coaplenent the tiae-dep^dent aodel technique and 
offer significant economic advantages when used in a judicious 
Banner. It is envisioned that the two techniques would be used 
in the following manner. Clean surge line determinations would 
be accomplished using the linearized stability prediction tech- 
nique. Frequencies at which amplifications occurred in a coa- 
pressor coaponoit/ducting system would be determined by conducting 
linearized frequoicy response analyses over the frequen~y range 
of interest. At points where amplification is shown to occur, 
the time-dependent aodel would be employed with large amplitude 
boundary conditions to determine harmonic content, waveform shape, 
compression component planar-wave- transfer coefficients, etc. 

Nhe-e a two sector parallel coapressor can provide an adequate 
estimate of loss in surge pressure ratio due to circumferential 
distortions, the linearized analysis stability prediction technique 
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would be enployed. If It was dooacd necessary to use a oulti- 
sectored parallel coapressor aodel, then it would be necessary to 
evaluate whether *’>*e linearised st^ility prediction technique or 
the tia»-dependent aodel offered econoaic advantages and to wake 
the choice for usage accordingly. 

The results of this study suggested areas that should receive further 
study. It is recosaanded that investigations be iapleaented to deteraine: 
1) The validity of the linearised analysis stability aodel and/or its 
dependence on the nuaber of stages operating on the positive portion of 
their characteristics, and 2) The potential existence of a physical inter- 
pretation of the eigenvalues that would be of benefit to tia coapressor 
designer. 
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AFPEMDIX A 


JADOBMB MAT&IX yOBHDIATIOW 


It was shoHD that the Jacobian (jh) Matrix (Equation 27) is a systen of 
first partial deriwatiwes of the ims» no ncntun and energy nacrobalances with 
respect to the state variables H, and s at a steady-state operating point. 
Rxani nation of the fomulation of these partial derivatives revealed that 
they could be expressed explicity as a function of the state varidiles, thus, 
elininating the need for finite difference approxinatioos. This elininates 
the dependence of the solution on perturbation size and therd>y, elininates 
a potential source of error. 


EafJi station requires at least three knoun flow properties to define the 
state of the fluid. At an internal station, the station values are obtained 
by interpolating between adjacent volune averaged properties. The fc-th 
volune fluid state is represented by the nacrobalances of nnnn. nonentim aid 
energy (Equations 1, 2 and 3). A typical volune (k) vlfh anropriate station 
notation (i) is shown in the following sdMaatic: 


I 


i-1 1 

Vl 

i-i . 

Vi 

1-1 1 

*k-l 


I 

I 

I 

I 




•* 1+1 

*1+1 

® 1+1 


I 

I 

I 


The set of three nacrobalances for eadi voliaa is 
and f^, respectively, as shown by Equation 25. 


denoted by 
In these equations: 




(pj.Hj) 


(41) 




(42) 


For stators. 




^ i-1 **1-1 *** i **1 i+1 **i+l^ 


(43) 


For rotors. 





(44) 
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For free volunes. 


r_ =■ 0 (45) 

®k 


For rotors. 

For free voloKS and stators, 

Sy = 0 (47) 


The A aatrlx is of order 3k for a fc voluae coopressor aodel. Ihe partial 
derivatives of the functions f-. are tatoi with respect to the voluae-averaged 
state variables contained in ilS equation. Figure 33 shows a schenatic repie~ 
sentation of Vol«nes 4, 5 and 6 of the J85-13 conpressor aodel. Ttese volimes 
represent the IGV-trailing-edge free volune, stage 1 rotor and stage 1 stator, 
respectively. The continuity equation, in each case, has derivatives taken 
with respect to and '^l+l. 



(48) 



aw 


i+i 



(49) 


The BoaenCUB equaticm provides derivatives of all surroundiog state 
variables ^ * 


^ !!a,. . 

»Pi “ ■ ^Ja^ apj ■ apj ^^i'Vi^ ap^ 


(50) 




( 51 ) 
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Continuity 


■i >0 *> *«it *^31 Ci *3k «| • * 
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( 52 ) 


3W. 3P. 3F_ 

t O i A 5 f A -A \ s 

38^ 38^ i 3Sj i+l' 3Sj 


!!k 

3p 


i+1 


3P 


i+1 


" (A.-A_,) + ® 


i+1 ‘*i+l^i+l*o *‘*i+l ^‘’l+l * ***■ ^i+1 


(53) 


»“k 


2W. 


i+1 


3F, 


B 


3U 


i+1 


‘*i+l^i+l®o *“i+l 


(54) 


3W^ 


DP 


i+1 


3P^ 


3F, 


38 


(A,-A,^,) + 


i+1 


»«i+l *«i+l ‘ **"1 »«i+l 


(55) 


Two additional derivatives are needed for stator volunes. 


®pi-i ^‘’i-i 


(56) 


3W. 


3W 


i-1 


3F. 


3W 


i-1 


(57) 


The energy equation provides derivatives taken with respect to s. and 
However, for rotor voluaies, the entropy production tern (Sp) yields 
additional derivatives with respect to Px>^i> ^1+1* 

For free voluses and stators, (sx<fX ~ ^i^ “ 0 




3s 

3s 


k 

i+1 



(58) 


(59) 


66 



For rotor voluses; ~ s^) 0 


>»! ■ 4 *^ [■ 


( 60 ) 


(61) 


3Si P,V^ [ 1 3s J 








(62) 


(63) 


3s. 


*®i+l »»i\ 


(64) 


A general equation expressing static pressure in terms of static density 
and entropy can be easily derived. 


P - P , ( ^ 
ref I p 


(t)’ 


e (s - s^^^Ht-D/R 


(65) 


Derivatives of Equation 65 with respect to density and entropy are: 


3p “ p 

(66) 

3P P(y- 1) 
3s * R 

(67) 
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formation of the derivative of 


( 68 ) 

( 69 ) 

( 70 ) 

( 71 ) 

( 72 ) 

( 73 ) 

( 74 ) 

( 75 ) 


parameters associated with the 



Equations 5 throu^ 10 are expressions used for decemlning the blade 
force tern (Pb)> An expression relating absolute flow angle (a) to relative 
flow angle (6) In terms of local state variables can be derived. 


tan 0 « - tan a (76) 

W 

Derivatives of the above equation and other fundamental relationships 
have been combined in the formulation of the blade force term derivatives 
and are shown In their final form below. 

For stators: 



p a (tan Bca) 

^i" 



- a (tan B«) 

awj " 

T awj 



B 

n tan 1 

®pi+l 

’“i+i ■ 




+ T — tan Bea 
3Pi 


awi 






3Pl_l ^ 3Pi-i ®Pi-i 


Boo 


.. a (tan 6») . 

_ = + _ tan 


6<» 


For rotors: 


3F, 


FpT-^T 


a (tan 6.) . T 
V. P»/ . + ^ tan Boa - 


3F, 


DZ 


3p 


3P. 


ap. 


(77) 

(78) 

(79) 

(80) 
(81) 

(82) 

(83) 
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(84) 


!!> . F + !!i 


!h 

de. 


3F, 


DZ 


3s. 


(85) 


!!b_'!!i_ 

^‘*1+1 ° ®Pi+i 


tan Bo, 


( 86 ) 


3F. 


3F„ 




tan 3« 


(87) 


Equation 12 is an expression for the entropy production term (Sp) in 
the energy equation and represents the blade row losses in rotor volumes. 



( 88 ) 


(89) 


(90) 


(91) 


( 92 ) 
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Hence, Equations 48 through 93 completely define the elements of the Jacobian 
matrix. 

Boun**ary conditions were Imposed at the ccxnpressor-model inlet and exit 
volumes. It was previously mentioned that each station requires at least 
three known flow properties to define the state of the fluid. Model-inlet 
total pressure and total temperature are inq>osed explicitly as boundary con- 
ditions, while the density is lnq>lled. At the model exit flow functfon was 
specified in a manner analogous to specifying a choked exit boundary condition, 
while the flow and entropy were implied. 

A parallel compressor model based upon an extension of the Dynamic Blade 
Row Compression Componert Stability Model was used for the dlstorted-inlet- 
flow studies. Parallel compressor theory is based on the assumption that 
each sector operates on its cloan inlet characteristics and exits to a uni- 
form static-pressure interface downstream of the conqtressor exit guide vanes. 
Flow from each sector is summed prior to entering the downstream volume and 
the entropy is implied. Inlet and exit boundary conditions are identical to 
those used in the clean- Inlet-flow model. 
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APPENDIX B 


frequency response input matrix formulation 


The analytical formulation of the partial derivatives of the input 
response (B) matrix is analogous to the methods, previously described in 
Appendix A. For a k - volume system, partial derivatives of the macrobalances 
of mass, mcxnentum and energy are taken with respect to each input (u„) while 
holding all state variables constant. The resulting matrix xs of the form 







'^3k 

3u^ 
L_ ^ 

3u 

n 


(94) 


The two Inputs used in this study are the compressor inlet total pressure 
and the compressor exit flow function. C''.isequently , the inlet and exit free 
volumes, respectively, are the only voi .mes containing these input parameters. 
All ^ matrix elements associated with the remaining volumes are zero. 

The first column of the ^ matrix contains the partial derivatives of the 
state equations with respect to inlet total pressure. The continuity and 
momentum equations (Equations 1 and 2), of the first volume, provide the only 
non-zero elements for the input total-pressure oscillation. In these 
equations for this particular volume, station subscript "1" corresponds to the 
compressor Inlet station, where 


'^i ' ”inlet ^^tinlet^ 

(95) 

^i ° ^inlet ^^*^inlet^ 

(96) 

Expressions for the ab !ve equations can be easily derived, 
partial derivatives of the static equations are 

The resulting 

• 

9pl j 

3Pti 

(‘»7) 
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dw 

3P 


ti* St [kw 




3P. 


^ “ 3Ptj " ^i+1^ 


) 


(98) 


'is. 


0 


(99) 


The second coluan of the ^ Batriz contain.? derivatives with respect to 
exit flow function. The BoaentuB equation asst.'ciated with the last free 
voluBe provides the only non-zero eleaents. At the teralnating station, 
static density and pressure are expressed in tens of the local flow function. 


P 

P 


i+1 ^exit 

(FT) 

(100) 

i+l ^extt 

(FF) 

(101) 


All aatrix eleaents associated with the cootinolty and enexg 7 equations 
are zero. The resulting partial derivatives are of the fore 


3FF 


• 0 


(102) 


^=!or 


i+1 


do 


3P 

i+1 i+1 


dFF L. I 2 . 

^^^0**!+! ^i+1 


dFF dFF i+1 


dFF ^*i ■ ^i+l^J 


(103) 


!!i 

3FF 


0 


(lOA) 


Once the input response aatrix is generated, all scale factors applied to 
the A aacrix ausr be applied to the B aatrix and the frequencies used to 
evaluate toe fn^iuency response characteristics. 
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APPBIDIX C 


DESCRIPTIOW OF MATRIX TRAMSFER FUHCTIOM SOLOTIOII TECmiqPE 


Given the A, and C natrices (Equations 28 and 29) of appropriate 
diaension, the problea is to obtain the aatriz transfer function given by 
Equation 33 

- C(SI - A)-l . 

The technique that has been prograaaed in Reference 8 considers only one 
eleaent of the aatrix at a tiae. For exaaple, consider the scalar transfer 
function between the i-th eleaeut of R (r) and the J-th eleaent of U (u). 


- = C. (si - A)"^ B. 
u — 1 — — — 1 


(105) 


where ^ is the i-th row of C and is the j-th coluan of B. 

The problea is siaplified by using the transfora given in Reference 16 
since it %rill eliainate the need for a Ci vector by asking tlM output r a 
state variable in the transforaed systea. 


I = I ®j (106) 

G = T Y (107) 

E = T A (108) 

such that 

G •• JE G + F u (109) 


where T is the identify aatrix except that the MC-th row is overwritten by the 
vector. The integer HC is the position of the eleaent of with the 
maxiaun absolute va.'ue. The aatrix is siailar to Inverses encountered in 
the proof of the Danilevsky aethod (see Reference 9) and can be written down 
explicitly. The output r of the unaodified systea is the HC-th SKxlified state 
variable Itius Equation 103 b coaes: 
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“ “ ( a I - E)‘^ F (110) 

or 

(s ^ - ^)r = F u (111) 


Thia is a aatrix of linear equations and in theory can be solved by 
eraser's rule. The solution is 


IT ^ det (W) 

u det (s ^ - ^) 


( 112 ) 


.There W is the systea Uronskian obtained by replacing the HC-th colunn of 
(s ^-E) with the F vector. The deterainant of (s ^-^) equals the deterai- 
nant of (s jl-A) and is the characteristic polynoaial of the sqrstea A aatrix. 
The Danilevsky aethod provides an efficient algoritha for obtainieg the 
characteristic polynoaial (see Reference 9). It uses siailarity tmsfor- 
nations to put the sutrix in a phase variable foza where a row of the 
cransforaed A aatrix becoaes the chr acterlstlc polynoalnal. The deterai- 
nant of W is found using a clewer algoritha developed by Eollinger (Refer- 
ence 17). Using a few properties of deterainants. be derived a aethod for 
evaluating the deterainant of tf using two evaluations of the Danilevsky 
type expressions: (1) The denoainator polynoaial is the characteristic 
equation of the systea transfer function and is obtained by the Danilevsky 
aethod using the A or E aatrix and (2) The nuaerator polynoaial is obtained 
in two o»re characteristic polynoaial calculations using Dollinger's aethod. 
The first characteristic polynoaial is obtained by v^lacing the HC-th 
col«m of £ by -F to obtain Z2; the second characteristic polynoaial is 
obtained through use of a aartix of order (n-1) and is obtained by deleting 
the MC-th row and coliam froa E to obtain Z3. The oxaerator polynoaial is 
coaputed as 


det (W) « Z2 - Z3 (113) 

Thus, systea transfer function is writtoi as 


o (W) 

det (s 1 - A) 
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